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, The problem of dark energy arises due to its self-gravitating properties. Therefore explaining 

' vacuum energy may become a question for the realm of quantum gravity, that can be addressed 

CNJ I within string theory context. In this talk I concentrate on a recent, string-inspired model, that 

^ ' relies on nonlinear physics of short-distance perturbation modes, for explaining dark energy 

O I without any fine-tuning. Dark energy can be observationally probed by its equation of state, 

' w. Different models predict different types of equations of state and string-inspired ones have 

\ a time dependent w(z) as their unique signature. Exploring the link between dark energy and 
string theory may provide indirect evidence for the latter, by means of precision cosmology data. 



^ '. 1 Introduction 
(N 



(N 
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There is still no fundamental physical theory of the very early universe which addresses 
issues that arise from the regime of transplanckian physics. It is the lack of a fundamental 
theory, valid at all energies, that makes the model building of the transplanckian regime 
^ ■ very interesting. The main issue is how much are the known observables affected by the 
unknown theory. The apparently ad hoc modification of the dispersion relation at high 
Qh! energies is contrained by the criterion that its low energy predictions do no conflict the 
^ I observables. 

In this talk we address the possibility that transplanckian regime can contribute to the 
' observed dark energy of the universe through a counterintuitive UV/IR mixing|l[]. The 
I physics mechanism that gives rise to dark energy is the freeze-out of ultralow frequency 
short distance modes by the expansion of the background universe. This transplanckian 
dark energy (TDE) model is motivated from closed string theory in a toroidal geometry by 
invoking superstring duality. 

The transition from string theory to conventional cosmology is becoming increasingly 
important to theoretical physics. As we describe in 0] corrections to short distance physics 
due to the nonlocal nature of strings contribute to dark energy. The possibility to detect 
their signature observationally is thus very intriguing. 

A potential way for detecting dark energy and distinguishing it from a pure vacuum 
energy A is through its equation of state. I will briefly touch upon the calculation of the 
stress-energy tensor, in order to address the puzzle of "cosmic coincidence" (see for 
details), and the time dependence in the equation of state w{z) predicted by this model. 



2 The TDE Model 



Let us start with the generahzed Friedmann-Lemaitre-Robertson- Walker (FLRW) hne- 
element in the presence of scalar and tensor perturbations. For simplicity, we will take 
the class of inflationary scenarios that has a power law solution for the scale factor a{7]), 
a{ri) = \ric/i]\^, where /? > 1 and Ir/d = l^/Hij]^. The initial power spectrum of the per- 
turbations can be computed once we solve the time-dependent equations in the scalar and 
tensor sector. The mode equations for both sectors reduce [|, || to a Klein-Gordon equation 
of the form 



a" 



/i„ = . 



where the prime denotes derivative with respect to conformal time, and n is the comoving 
wavenumber. Therefore, studying perturbations in a FLRW background is equivalent to 
solving the mode equations for a scalar field related (through Bardeen variables) to the 
perturbation field in the expanding background. The above equation represents a linear 
dispersion relation for the frequency a;. 



-^=P^ = ^- (2) 
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The dispersion relation of Eq. @) holds for values of momentum smaller than the Planck 
scale Mp. There is no reason to believe that it remains linear at ultra-high energies larger 
than Mp. 

In what follows, we replace the linear relation oj'^{p) = = r? jaij])^ with a nonlinear 
dispersion relation = F{p). Therefore, in Eq. (|l]), should be replaced by: 

nlff = aivrFipy = aivrF[n/aiv)]\ (3) 

Let us consider the following Epstein function |^ for the dispersion relation, in the 
squared bracket in Eq. (|l|): 

= np)=p'{T^ + ^ + 7T^^.\^ (4) 
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<ff - - = a\r^)F\n, v)=n^-^ + i^ + — ^ ) , (5) 



where x = {p/pcY^^ = A\t]\, with A = (I / \r]c\) (n/ pcY^ ^ , in general for power-law inflation 
with f3 > 1. This is the most general expression for this family of functions. For our 
purposes, we will constrain some of the parameters of the Epstein family in order to satisfy 
the features required for the dispersion relation as follows. 

Imposing the requirement of superstring duality fixes €2, and the condition of a nearly 
linear dispersion relation for p < pc requires that 

^2 = 0, 1 + 1 = 1- (6) 
Still we will have a whole family of functions parametrised by the constant ei, as can be 




seen in Fig. |l]. 

The initial vacuum state is well-defined in this model and it is a Bunch-Davis vacuum, 
due to the dispersion function going asymptotically flat at early times. At late times the 
solution becomes a squeezed state by mixing of positive and negative frequencies: 

H'n ri—*+oo i ~r , o , j 



with |/?nP being the Bogoliubov coefficient equal to the particle creation number per mode 
n, and ~ ^Jeln. Using the linear transformation properties of hypergeometric functions 
we obtain: 

1 smhzTT^/ei 



Clearly the spectrum of created particles is nearly thermal to high accuracy. 

We define the tail as the range of those transplanckian modes in fig.l whose frequency 
is less or at most equal to the present Hubble rate, (see Fig. 2). Since H has been 
a decreasing function of time, many modes, even those in the ultralow frequency range, 
have become dynamic and redshifted away one by one, everytime the above condition is 
broken, i.e. when the expansion rate H dropped below their frequency. Clearly, what is 
left from the tail modes at the present Hubble scale have always been frozen. They contain 
vacuum energy of very short distance, hence of very low energy. The last mode in the tail 
(at infinite momenta) would decay when and only if if = by acquiring a kinetic term. 
The tail starts from some value "pu which must be found by solving the equation 



(9) 
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Fi gure 2: The range of modes in the tail, pn < p < oo, defined by Eq. Hq is the present value of the 
Hubble constant. 

The range of the modes defining the tail is then for pn < p < oo. Their equation of 
state depends on the evolution of H and is a comphcated tracking solution because they 
contribute to the expansion rate for H and their equations of motion are coupled to the 
Friedmann equation for expansion. 

However let us first calculate their energy contribution today, by asking which "tail" 
modes have always been frozen from infiationary times to present and how much energy do 
they contain. The energy for the tail is given by: 

{Ptau)=^-^ r^Uj{p)p'dp, (10) 
277^ JpH P 

while the expression for the total energy is: 

(Pi0tal) = ^j^ —^Uj{p)pdp. (11) 

The numerical calculation of the tail energy produced the following result: for random 
different values of the free parameters, the dark energy of the tail is ptaii = 10~^^^/(ei), 
times less than the total energy during inflation, i.e. -^^^ = 10~^^^/(ei) at Planck time. 

Ptotal 

The prefactor /(ei), which depends weakly on the parameter of the dispersion family ei, 
is a small number between 1 to 9, which clearly can contribute at the most by 1 order 
of magnitude. This is an amazing result! The lack of fine-tuning in obtaining the right 
ammount of dark energy can be understood from the following feature: due to the decaying 
exponential, the main contribution to the energy integral in Eq. ([ToD comes from the highest 
value of this exponentially decaying frequency, which is the value of the integrand at the 
tail starting point, pn ~ 0(Mp), i.e., 

^j^^-w/(^"^-wr''''^- '''' 



From the physical requirement that the tail modes must have always been frozen, the tail 
starting frequency uj{ph) is then proportional to the current value of Hubble rate Hq (Eq. 



3 TDE Dispersion Relation from Closed Strings in Toroidal Cos- 
mology 

Here we attempt to motivate the exponentially supressed transplanckian dispersion function 
from T-duality of the following string theory model. Let us consider the Brandenberger-Vafa 
model P] of a D-dimensional anisotropic torus with radius Ri, by including the dynamics 
of both modes: momentum modes, pi^i = m/Ri (where m is the wavenumber), and winding 
modes with momenta p2,i = wRi/a'. The dimensionless quantity for the radius is Ri = 
Rij^^ ^ where a' is the string scale. Based on the arguments presented in p[, we choose a 
cosmology with three toroidal radii equal and large i? 3> 1 in units of the string or Planckian 
scale, with the other {D — 3) toroidal radii equal and small Rc <^ 1. Here the subscript C 
refers to compactified dimensions. Then, R{t) becomes the scale factor for the 3+1 metric 
in conventional FRW cosmology R{t) = a{t), while Rc corresponds to the radius, in this 
factorizable metric, of the D — 3 compact dimensions zj that decrease with time, 

dsl = -dt^ + 4nR^{t)dx'^i + ^'KRl{t)dz] = a{rif[-d7f + dy'^] + ds]j_^ . (13) 

These solutions were found by Mueller, Ref. 0, by ignoring the backreaction of string 
matter on the geometry. 

The partition function for this system was calculated, from first principles, by summing 



up over their momenta in 



Z = ^e-"'^^% (14) 
where n^- is the number of strings in state a with energy e^- 



= Po = y (^) ' + (wRy + N + N-2, (15) 

and a counts over (m, w), with the constraint N — N = mw for closed strings where and 
are the sums over the left- and right- mover string excitations, respectively. By now, in 
Eq.(|r^), we are considering only the large 3 spatial dimensions. The string state can also 
be described by its left and right momenta, ki = pi +p2, = pi —p2- The string state for 
left and right modes can be expanded in terms of the creation and annihilation operators 
am, On, with higher excitation string states given by iV = X^^i Q^-nttn (similarly for N), 
and string energy Lq + Lq = p\ + p\ + {N + N — 2) / a' . 

We would like to write the path integral for this configuration in terms of quantum 
fieldsQ The path integral is calculated from the hamiltonian density. The hamiltonian 
density over the fields in configuration space is extracted from the string spectrum in such 

^Below we use quantum string equations under the assumption that the dilaton is massive and stable. 



a way that its Fourier transform in fc-space corresponds to the string energy expression Eq. 

m- 

Our system of winding and momentum modes is described by nonequihbrium dynamics 
due to the expanding background spacetime. All the information about the evolution of 
these modes will be contained in the effective action. The kinetic terms are unambiguous 
while for the interaction terms we must appeal to simplicity and the requirement of T- 
duahty. 

Correlation functions are obtained by using the correspondence between the Euclidean 
path integral of the persistence vacuum amplitude \ {in\out)\'^ and the partition function Z. 
All the string quantum operators below are promoted to quantum field operators with the 
corresponding hamiltonian density 7i{t, x) in configuration space derived from the quantum 
string hamiltonian T-C{t). 

The following calculations are done in the conformally flat background through the 
scaling of the fields and operators with the conformal factor a{r]). The momentum field 
(f)i{R,x) and the winding field (j)2{R,x) are defined by the relation: 

= J e'P^^<j)Md^p^ , J \V(j)^\^d^x = J (fpdUPi) , (16) 

where 

W = Rd/dx =d/dy, (17) 

and Pi = pi, p2- Let us also define two new fields, iPl{R,x) and tpR^R^x), with momenta 
kii kn that are the left and right combinations of the Kaluza Klein momentum and winding 
modes 

^^l{R,x) = MR,x) + MR.x), (18) 
MR,x) = MR,x)-MR,x). (19) 

These fields live in the expanding (3+1) spacetime dimensions. 

The Hamiltonian density ansatz that would describe the energy of our two string states 
in the D = 3 expanding dimensions with energy H = Lo + Lq, including the oscillators 
from string's higher excitations (A^ + N — 2)/a', is similar to the hamiltonian of spin waves 
in a periodic lattice^. Our lattice spacing is given by the string scale y/a'. Therefore the 
hamiltonian density can be written for this dual lattice in terms of wave functional "spin" 
fields ^l{R,x),^r{R,x) of Eqs. (|l|), (|l|) as follows 

(20) 

where the fields iIjl, ipR are expanded in terms of the mode functions u„,Mn, 
^Torus is obtained by identifying the first and the last lattice sites, thus the periodicity. 



and bn, bn are the normalised quantum creation and annihilation operators of a„, (in- The 
commutation relation for the unnormalised operators are such that [«„, a+] = uj±6nm with 
uj± the frequency of left and right moving modes. 

The periodic lattice condition N — N = mw introduces an interaction term in the 
hamiltonian H3 of the form V^'lVV'/j. In terms of the 2-component state "^en = {'ipL, 'ipR), 
the hamiltonian reads, 

Hs = \V'^EN\^ + V<i/ENXV'^EN+ml\<i/EN\^ + 9l\'^EN\^ + i92-2gi)\^ENX<i/ENW (22) 

with 

The system is known as the dual momentum-space lattice, and for g2 = 2gi reduces to 
the XYZ model of condensed matter. Let us for simplicity limit to the XYZ model case, 
92 = 2(71, for the rest of this talk. 

These periodic lattice systems studied in 3+1 dimensions in terms of Bloch wavefunc- 
tions have a solution which respect lattice translation invariance, exp{—pl), with the lattice 
spacing "/" equal to the string scale \/^. The interaction term, in the tight-binding approx- 
imation, lifts the degeneracy between the energy eigenstates due to the leakage/tunnelling 
of the wavefunction from one lattice site to the neighbour site. As a result the gap energy 
produced between the ground (bound) state and higher excitation states is 

p^Ap = / I cos(2^)| = p2 \2cos'^{e) - 1| , (24) 

in which 



pl = pV^= ^a'ipl+pl) = sjC^Y + (wRy , (25) 
and 6 ^ 6 + ipl. Therefore, 

Ap < 2 cosh^ (pi) -1. (26) 

The first term in 6' is a pure phase of rotation of the "spin-wave" in the dual lattice, but 
the second term describes the tunnelling of the wavefunction to the nearest neighbour^. 
The gap energy of Eq. (24) introduces a correction to the kinetic energy, such that in 
momentum space the hamiltonian reads 

Hs = ZpP^l'^ENl^ + mll'^ENl^ + 9i\'^en\^ , (27) 

with Zp = 1 + Ap. One can evaluate the reduced 3 dimensional partition function Z from 
the above hamiltonian [Q. We want to separate our modes into system (S) + environment 
(E) degrees of freedom, and coarse grain by integrating out the degrees of freedom for 
the environment. This amounts to finding out the backreaction of the coarse grained 
environment on the system, and eventually leads to the RGE'sQ. We will consider as 
environment all the short wavelength modes with momenta 

(E): T<P'' = -^IMR)^ + {wR)^Y^^ < A , (28) 
b J a' 



^ In condensed matter this is known as Coulomb dipole type of vortex interaction. 
^The running of the couphng constants with time depends on how one selects the environment and the 
system. 



where the cutoff A = [a')'^^"^ is the string scale because {a'Y^'^ is identified with the 
lattice spacing /, and h = a(t)/a{tQ) is the coarse grain scaling parameter, where to is the 
initial time. The scale factor a{t) plays the role of the collective coordinate describing the 
environmental degrees of freedom. Time in this procedure is playing the role of a scaling 
parameter and dynamics is being replaced by scaling, an artificial procedure, known as 
Kadanoff-Migdal transform, that relates the microscopic and macroscopic properties of a 
system based on the existence of scaling properties of the system in the infrared limit. 

The system modes are the ones with: 

(S): / < ^ , (29) 



^From the above definitions of system and environment, Eqs. (29) and (^), at initial times 
when ^ 1 we have p < A/b thus all our modes, momentum and winding, are in the system; 
but at later times when b > 1 more and more winding modes systematically transfer to 
the environment because the condition of Eq. (|9D, p = -^[{m/Rf + {wRfY'^ < A/b 
is satisfied only for vanishingly small winding numbers w ^ 0. As t becomes large, the 
system contains m < RA, w = 0, i.e. all the modes except m < RA,w = have been 
transfered to the environment. 

After integrating out the high energy environmental modes in the above action, one 
can extrapolate the correlation function from the coarse grained effective action of the 
system. The canonical two-point correlation function at high energy for system-environment 
interaction is calculated from the path integral of the canonical fields '^s,e in momentum 
space (Fourier transform of G^^^). It is related to the correlation function of the original 
fields (which decreases at high energy) as follows 

{^e^e) = . (30) 

Zp 

where Zp = 1 + Ap and Ap is given in Eq. (pB]) . This is the crucial result for the interpre- 
tation of the cosmological dark energy. Because of the mass gap, the correlation function 
is suppressed exponentially in p-space. It is very familiar that a mass gap leads to an 
exponential fall off in x-space, but here for the dual lattice the exponential fall off is in 
momentum space. This may be traced to the T-duality of the closed strings and the resul- 
tant interchange of the IR/UV limits. Thus the two point function of the original fields (by 
dividing with the Zp normalisations factor), goes as an exponential for large momentum p, 
and polynomial for low momenta. 

The observed small value A ~ 10~^^° in natural units has an explanation in the toroidal 
cosmology of closed strings and thus the dark energy provides an exciting opportunity to 
connect string theory to precision cosmology. We may argue that numerically the size of the 
cosmological constant in the present approach is a combination of the string scale and the 
Hubble expansion rate in the sense that A/Mp^^^^^ ^ 10^^^° ~ (Ho/Mpianck)^- Therefore 
the correct amount of dark energy obtained by this frequency dispersion function does not 
require any fine tuning and relies, besides a physical mechanism (such as freeze-out), only 
on the string scale as the parameter of the theory. 



4 TDE Equation of State 



Nonlinearity of short distance physics and the breaking of Lorentz invariance results in a 
violation of Bianchi identity for transplanckian models. Therefore one needs to modify the 
stress-energy tensor [T^^,) in the Einstein equations, such that the modified ones satisfy 
Bianchi identity. 

Based on the equation of motion as our sole information for short-distance physics, we 
therefore use a kinetic theory approach for modifying Einstein equations in the absence 
of an effective lagrangian description. The assumption made is that a kinetic theory de- 
scription of the cosmological fluid is valid even in the transplanckian regime. Despite its 
nonlinear behavior at short distances, this imperfect fluid shares the same symmetries, 
namely homogeneity and isotropy, as the background FRW universe. Then the corrections 
T^i, to the stress energy tensor T^,^ will also be of a diagonal form 

V = (e + n)M^u^ + Ug^^ (31) 



In a similar manner to particle creation cases in imperfect fluids [|T2[, the highly nontrivial 
time-dependence of the mode pn and the transfer of energy between regions, due to the 
defrosting of this mode across the boundary p^, gives rise to pressure corrections in the fluid 
energy conservation law. The defrosting of the modes results in a time-dependent "particle 
number" for regions near pn- From kinetic theory we know that this "particle creation", 
(the defrosting of the modes), gives rise to effective viscous pressure modifications fl^. The 



term 11 denotes the effective viscous pressure modification to the "bare" pressure, {p). 

The criteria we will use for modifying T^,y is that Bianchi identity must be satisfied [|13 
with the new expressions for pressure, P, 

E[p + 3H{p + p + n)] = S,[p + 3H{p + P)] = , (32) 

The presence of pressure corrections, Ilj, (where i = II, H denotes transplanckian defrosted 
and tail modes respectively) is due to the exchange of energy between the two regions, 
from the defrosting of the modes pn at the boundary. This is directly related to the 
time dependence of the boundary ph, which in turn is going to be controlled by the Hubble 
parameter H. In essence, there is an exchange of modes between the two regions. Although 
the number of particles^ in each of these regions, Njj and N^, is not conserved, their rate 
of change, in the physical FRW Universe, is related through the conservation of the total 
number of particles which contains both of these components 

Nt = 0, (33) 



The contribution terms to pressure, Ilj, are related to decay rates through |[T2[| 



3HUh = ~{ph+Ph)Th, (34) 



^We are loosely using the term particle here to refer to the wavepackets of the transplanckian modes, 
centered around a momenta pi. 



Details of the calculation can be found inP|. The value of the tail's decay rate Th is: 



TH^3H+^^3Hi^-^^), (35) 

where Wtotai = Ptotai/ Ptotai- Notice that Th is positive for all equations of state Wtotai < 1 
and thus it slows down the dilution of the tail with the scale factor. 

The effective equation of state resulting from the modification Th of stress-energy tensor 
becomes: 

WH = (!^^^^^|^) . (36) 

This is a tracking solution. During a radiation dominated universe, wh = — 1/3 and only 
after matter domination, the tail starts acquiring a strongly negative equation of state 
until it asymptotically aproaches wh = — 1 when it becomes the only component left in the 
expansion equation. Therefore cosmic coincidence is explained naturally by the intrinsic 
time evolution of the tail which tracks the background equation of state Wtotai through its 
coupling to the Hubble parameter. 

Careful analysis of the highly anticipated new data from CMB, LSS and SNla will reveal 
more scrutinity between various mo dels [p!^ . 
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